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Abstract 



The problem of dark energy(DE) is the greatest challenge for modern physics, and 
therefore this thesis is dedicated to the modeling establishment of DE. In the first 
chapter, qualitative and descriptive methods are employed to draw a rough map for the 
cosmos in its finely analyzed era. Then the dynamical equations of cosmology based 
on the general relativity are introduced, together with the measuring methods and 
parameters. The second chapter reviews the results of cosmic microwave background 
radiation and supernova observations, which conclude the fiat topology and accelerat- 
ing expansion of cosmos, and give rise to existence of DE and the associated dynamical 
mechanism problem. Taking WMAP data as an example, we introduced the practical 
parameters. In the end we classify the existing DE models into three categories for 
further discussion. The first two chapters make up the basis for the whole thesis. 

Chapter 3 reviews the fundamental structures and research progress of two scalar 
fields, Quintessence and Phantom, which correspond to the case of EOS beyond -1 
and below -1 respectively. Chapter 4 reviews the fundamental structures and research 
progress of another two scalar fields. Phantom and Quintom with non-canonical la- 
grangian, both of which have EOS crossing -1. But Phantom field is real, while Quin- 
tom complex. Chapter 5 reviews two non-scalar-field DE formalism of K-essence and 
Chaplygin gas, which are both rooted in modern quantum field theory. Totally six kinds 
of DE is analyzed, all of which have positive energy density, negative pressure and could 
drive the lately accelerating expansion. And they are compared and identified in details. 

New work gathers in Chapter 6 and Chapter 7. Firstly we make a comprehensive 
analysis of the existing work, pay particular attention to Gondolo and Preese's idea of 
treating Cardassian energy term as relativistic perfect fiuid(GF fluid), point out that 
a potential Cardassian term should meet three conditions, and review three existing 
Cardassian terms (power style, polytropic style and its modification, exponential style 
and its modification), and eventually put forward the newly found hyperbolic cosecant 
Cardassian. 

Then the Cardassian dynamical cqTiations are introduced generally and logically 
under GF fluid scenario, together with the flowing process of constructing phase space 
and differential dynamical systems from Fricdniann equation. Hyperbolic cosecant Car- 
dassian term is employed for concrete computing. The analysis proceeds in two cases, 
namely a unifled description of matter and radiation energy density (case 1) and a sep- 
arate description of matter and radiation terms (case 2). Formalism of case 2 is more 
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exact at the expense of more complicatedness, and due to the mathematical symmetry 
of matter term and radiation term in hyperbolic cosecant function, the differential dy- 
namical equations are considerably simplified. Phase space and dynamical systems for 
both cases are achieved. When we calculate the critical points for case 2, amazingly 
interesting behaviors of self-consistency and auto-normalization are exhibited, which is 
a strong support for the new model, along with a forever positive sound speed. 

The process of virial collapse in Cardassian cosmos is analyzed. Power-style Car- 
dassian term is employed for its simplicity. Calculation declares that virial collapse of 
matter alone is forbidden. Yet Cardassian has excellent ability for virial collapse, after 
the virial collapse ending up with a stable sphere, the ratio of the ultimate radius to 
the original radius depends on the adjustable parameters in Cardassian term. And, the 
mixture of GF fluid and matter could conduct virial collapse, the ratio of the ultimate 
radius to the original radius depends on the adjustable parameters in Cardassian term, 
too. No singularity is generated. 

The creative work in this thesis incudes the introduction of hyperbolic cosecant 
Cardassian, which is the fourth ever available Cardassian term ever found, after the in- 
troduction of exponential style in 2005. This work, along with the calculating of power 
Cardassian virialization, enriches the research of Cardassian cosmology. 

Key Words: Friedmann Equation, Dark Energy, Scalar Field, Cardassian, Hy- 
perbolic Cosecant Function, Virial Collapse 
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Chapter 1 

Introductions 



1.1 Fundamentals of Cosmology 

Cosmology, which ranks one of the most eneouraging realms in natural science, is the 
scientific branch researching into the cosmos structure, origin and evolution on the base 
of astronomic data and physics principles. Modern cosmology is basically supported by 
three events, namely the cosmic principle, the expansion of the universe together with 
Hubble Law, and the Cosmic Microwave Background Radiation. 

1.1.1 the Cosmic Principle 

Shortly afterwards the establishment of General Relativity, A. Einstein computed 
the global behaviors of the cosmos. To stay in harmony with Mach Principle and to sim- 
plify the calculating process, he made a basic assumption, which is nowadays honored 
as the Cosmos Principle. It declares that the universe is homogeneous and isotropic 
on large scales. Under modern considerations, the so-called large scale means covering 
a distance of 10^ light year or more, which belongs to the cosmic scale. Homogeneous 
means different points identify with each other and could not be tell apart via physi- 
cal experiments, while isotropic refers to the fact that , when taking observations into 
arbitrary directions at any coordinates origin, the scenarios are the same. A most sig- 
nificant application of this principle is to simplify the form of the space-time metric 
quantitatively. At the present days, the Cosmic Principle has already been established 
firmly on the base of observation data. As a matter of a typical sample, the results 
of the Willkinson Microwave Anisotropic Probe (WMAP) indicate that the universe 
is in high homogeneity and isotropy, only with tiny fluctuations at the magnitude of 
10"^ ~ 10~^, which is in fact essentially required for the generation of galaxies according 
to the standard model of cosmology. 
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1.1.2 Cosmos Expansion &c Hubble 's Law 

In 1929, Hubble was engaged in comparing the light spectrum of far distant galaxies. 
To his surprise, he discovered a systematic redshift with the galaxies beyond our milky- 
way. As to redshift, it means the whole spectrum translates to the red wave interval. 
According to Doppler Effect, all galaxies are leaving us, and the spatial scale of the 
universe is undoubtedly expanding. After some algebra with the redshift amount and 
the distance to the galaxies, Hubble found the routing that the farther they are, the 
quicker they leaves, at a speed v proportional to the very distance r: 

V = Hor (1.1) 

where the coefficient Ho is honored Hubble Constant, which is simply the present-day 
value of the more generalized dynamical Hubble parameter. Currently, the results mainly 
lie in the interval Hq 60 ~ 70Km • • Mpc~^. 

1.1.3 Cosmic Microwave Background Radiation 

The background space of galaxies is a radiation field that is highly homogeneous 
and isotropic, whose power spectrum curve agrees with that of the 2.7K blackbody 
perfectly. The so-called Cosmic Microwave Background Radiation (CMBR) is a key 
tool to understand the structure, origin and evolutionary history of the universe. It was 
researching into CMBR that lead to the era of finely quantitative analysis of cosmology. 

CMBR coTild translate quite a lot of information on the young cosmos. For instance, 
the perfect homogeneity and isotropy suggest that the cosmos is also homogeneous and 
isotropic at its youth era, and the scattering surface is spherical, on which the gas are 
of identical temperature. The background photons are in thermal equilibrium before 
decoupling, thus the spectra, which describ how the radiation strength depending on 
the power at arbitrary direction, are all blackbody-like. The temperature fluctuations 
at present-day represent the medium density fluctuations at early time, which magnify 
via gravitation and form the clustering structure, or the inhomogeneity on the small 
scale. CMBR plays a pretty important role in the standard model of modern cosmology. 

1.1.4 Global Sceucirio 

The standard model of modern cosmology comprises two parts: One is the heat Big 
Bang theory, which was put forward by G. Gamow at around 1948, the other is the 
inflation theory, which is born in early 1980s when Guth,Linde, Albrecht and Steinhart 
employed this mechanism to beat the difficulties that Big Bang ran into at extremely 
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early era. 

The heat big bang theory is based on the observational evidences of Hubble Expan- 
sion and the Helium-4 amounts, and the theoretical formalism of General Relativity 
and solar heat nuclear reaction. It succeeds in explaining the genesis of light chemical 
elements, Hubble redshift and forecast the 5K CMBR (although the measured value 
being 2.7K, this prediction is still a greatest success of mankind intelligence). The 
big bang formalism accounts for the process from 10~2 second (T ~ lO^^K) to cur- 
rent time (~ 10^^,s,2.7K), but falls in trouble with the so-called singularity problem, 
the horizon problem, the asymmetry of baryonic numbers, the flat state problem and 
origin of galaxies, all of which are concerned with the extremely early era of the cosmos. 

The inflation theory is based on general relativity, as well as the standard model 
for high energy particle physics established at that time. It aims at the time scale from 
10^^ second dating back to 10~^^ second (Plank time), and solved all the problems the 
big bang theory come across all at once, except for the singularity problem, which is 
principally impossible to overcome under the framework of classically general relativity. 
There seems no avoiding the singularity problem, and only the quantum cosmology 
theory declares to possess the power to solve the diagnose. 

1.1.5 First three Minutesin Fine Cosmology 

The time scale is generated together with the Creation of Big Bang. When time 
reaches 10~^^ second, the scale grows to 10"'^^ cm, the four finds of elementary forces 
are unified as a grand force. Only energy lives, without any matter, and the world is 
dominated by quantum uncertainty principle. This era is called the Plank era or the 
Grand Unified erap]. 

when it goes to 10~'^^,s, temperature drops to 10^^ K, and certain kinds of particles 
are produced. Although some particles annihilate with their anti-particle companies, 
the creating rate is faster than the annihilate rate. Temperature slowing down, gravi- 
tation and the nuclear strong force decouple from the grand unified force one after one. 

Then since comes the inflation era. The spatial scale doubles ev- 

ery 10~^^ second. In a period of lO^^^s, the scale expanded by 10^'^'^ times, ending with 
a spatial scale of 3 x W~^^cm and temperature of 3 x W^^K. By the end of inflation, 
the electroweak force decouples into electromagnetic force and nuclear weak force. 
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When it comes to 10~ th second, the corresponding temperature is 10 if. The 
highest energy produced on the large hadron colhders is able to simulate this situation. 
And the validity of the numerical calculating and experimental simulating stands test- 
ing and win a wide agreement. 

When it comes to 10~^ second, the temperature decreases to W^^K, and "quark- 
gluon soup" is generated. Quarks and gluons are in free state originally. With the 
temperature continuously doing down, a gluon will bound three quarks together to 
form a nucleon (proton or neutron). 

When it comes to 1st second, temperature drops to 10 billion Kelvin. The con- 
tains of the cosmos are photons, quarks, electrons, neutrinos and other particles, and 
of course, proton and neutrons. Temperature decreasing, neutrons begin to decay into 
steady proton, ending with a proton- neutron number ratio 7:1. 

With time slipping down, the creation rate drops, and annihilation rate of particle- 
antiparticle pairs exceed that of creation. Large amounts of electron-positrons annihi- 
late to produce more photons. Eventually the annihilation era finished, the bulk being 
the remaining quarks and electrons (leptons). The cosmos evolves into a new creating 
era, and above all more protons and neutrons are created. 

At about 100 seconds after the big bang, it gets 1 billion Kelvin. Protons and 
neutrons can never escape from the constraints of nuclear strong force and forms Hf. 
Then turns to -ffe^, and iJe^ to -ffe|, with large quantities of neutrinos generated 
meanwhile. 

Three minutes after the big bang, the light nuclear synthesis process ends due to 
low pressure. Now about 25% of protons and neutrons forms Helium, with a few Hf, 
Lis, Be4, and the remaining neutrons decay into proton. This results in a universe of 
75%Hydrogen and 25% Helium. 

The standard model of particle physics also believes that, in the extremely early era, 
some other particles are produced, such as x-boson, Higgs boson. Weakly Interacting 
heavy particle, cosmic string and monomagneton .But none of these are explored. At 
three minutes after the big bang, the universe is filled with electron, nude hydrogen, 
nude helium and dark matter, as well as high energy 7. Because photons, which is the 
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medium particle of electromagnetic force, interact with free electrons and positrons via 
collision, absorption, emitting frequently, the lifetime and free-diatance of a photon is 
quite short, and information couldnot travel to the distance. So, the whole cosmos is 
in dark plasma state. 

0.3 million years after the big bang, it has reached a temperature as low as 4000K. 
Electrons are bounded in atoms, so photons could travel to distant places to compare 
notes. The cosmos is apparent now, and it expands steadily and silently, and enriching 
clustering structures are formed. The temperature of background radiation goes on 
falling ,to the present-day value of 2.7K. 

1.1.6 Dark Matter 

Just before the discovery of dark energy, the so-called dark matter is already a 
remarkable challenge to man's brains. Although Zwicky had suggested the existence 
of dark matter, it is not until the late 1970 did the idea is highly treasured. Now, it 
is commonly accepted the existence of dark matter. Lots of experiments lead to this 
result, the classical representative of which comes from the measured rotating speed 
curve of spinning galaxies. With the hard efforts of high energy physists, conclusions 
are drawn that the dark matter, which has hardly any interactions with photons, could- 
not be baryonic. Hence it is now generally called cold dark matter (hereafter CDM). 

There are several candidates of elementary particles for CDM. The most hopeful 
one is the weakly interacting massive particles (WIMP), whose mass is probably far 
beyond a proton. Among several WIMPs, a particular kind is the neutralino, with a 
typical mass of 102GeV, or 10~^^ gram. Other candidates include the axion with a 
typical mass of lO^^eV or 10^^^ gram, fuzzy CDM, neutrinos, etc[2]. 

Nevertheless, CDM might also be ancient black wholes that form before nuclearsyn- 
thesis, rather than elementary particles. And the lower mass limit of the blackholes is 
Hawking evaporating limit, namely 10^^ gram. 
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1.2 General Relativity & Cosmic Dynamics 



1.2.1 Gravitational Field Equation & FRW metric 

The research into cosmic behaviors on large scales is based on Einstein's gravita- 
tional field equation [1] [5]: 

G^u = R^^u - \gf.uR = SttGT^, (1.2) 

where natural units are employed to set c = /i = l,and G^i, is Einstein tesor, R^i, being 
Ricci tensor, i? being Ricci scalar, G being Newton's gravitational constant. 

Take another consideration, as for the four-dimensional spacetime that obeys the 
Cosmological Principle and is therefore homogeneous and isotropic, the metric is Friedmann- 
Robertson-Walker(FRW) metric: 

ds'^ = dt^ - a^{t) + ''^'^^^ + ^^"^ ^^'^^) (^-^^ 

where a{t) is the scale factor whose evolution is dependent on cosmic time i; r, 9, cj) is 
the set of comoving coordinates fixed on the medium, the medium expanding with the 
universe while the comoving coordinates remaining static, k refers to the curvature 
scalar. When A; > 0, the spatial topology is limited curved and namely closed universe; 
when A; = 0, the spatial topology is infinitely flat and namely flat universe; when k < 0, 
the spatial topology is infinitely curved and namely open universe. Moreover, if the 
unit for r is properly selected, the values of k can be reduced to he k = +1,0,-1 
correspondingly. Sometimes, FRW metric (II. 3p is transformed to take the formalism 
below for computing convenience: 

ds^ = dt" - a\t) [dx^ + fl{x){de^ + sin^ ed<i?)\ (1.4) 



where 



fkix) 



sin X, A; = +1 

X, k=0 (1.5) 



sinhx, k= -1 

1.2.2 Perfect Fluid Scenario & Conservation Laws 

Now let's combine the cosmological principle and the energy- momentum tensor 
which reflects the general states of matters, and get: 

T^u = ip + p)U^Uu + 'pg^,u (1-6) 
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where the energy density p and pressure p are dependent on t but independent of r, 6, 4>, 
and is the four-dimensional velocity vector: 



[/* = 1 (1.7) 

U' = (1.8) 

Nevertheless, the energy-momentum tensor for 4D perfect fluid also takes the form of 
Eq. (jl.6p . where p and p are measured in an inertial reference comoving with the fluid. 
Hence perfect fluid scenarios are generally employed as dynamical mechanism in cos- 
mology. 



The spatial media(taken as perfect fluid) are pretty homogeneous and isotropic, 
which determines the diagonal matrix of its energy-momentum tensor and the identifi- 
cation of the spatial elements. 

= diag{-p{t),p{t),p{t),p{t)), (1.9) 

The total energy-momenata ought to be conserving, hence there being the covariance 
equation: 

T^l"" = (1.10) 
as to = r, 9, <j) the above equation holds obviously; as for = O,one gets: 

d[a^{p + p)]= a^dp (1.11) 

and translate it into the formalism of the first law of thermodynamics: 

d{pa^) = -pd{a^) (1.12) 

which indicates, for a comoving volume element, the multiplication of its pressure and 
volume (i.e. volume work) pd{a'^) is the opposite value of the energy change of the mass 
element d{pa?). 

1.2.3 Priedmann Equation 

The expanding rate of the universe is reflected via Hubble parameter if, defined as: 

H = - (1.13) 

a 

Hubble parameter evolves according to observations, and is dynamically treated in the- 
ory too. The so-called Hubble Constant Hq refers to the presentday value of Hubble 
parameter in particular. 



12 



Under FRW meric,the non-null component for the Ricci tensor of gravitational field 
equation ()1.2p is: 

Roo = -3- (1.14) 

a 

R., = -(^ + 2^ + 2^)g,, (1.15) 



Ricci scalar being: 



a a 

the — O(time) component of the gravitation equation being 



This is the very Priedmann equation that represents the dynamical evolution of cosmic 
scalar factor, which indicates that that the expansion process depends on both energy 
term and geometrical curvature term. The i — i (spatial) component of gravitational 
equation (|1.2|) is: 

a a} k ^ , s 

2- + ^ + ^ = -SvrGp 1.18 

Meantime it is also induced to get Raychaudhuri equation which shows the evolution 
of Hubble parameter: 

ij = -47rG(p + p) + 4 (1-19) 

Eq ()1.12p . Eq (ll.l7p . Eq ()1.18p is connected via Bianchi identity, hence only two ones 
are independent. From p.lUp . (|1.17p . (|1.18p one could induce the energy- momenta 
conserving equation for cosmological perfect fluid: 

/9 + 3F(p + p) = (1.20) 

Unite (frT7l)(frT8]l to deduct \ and one would get a straight overlook on the cosmic 
expansion: 

a 47rG , , , , 

- = ^^ + 3^ 1.21 

Conclusion is directly drawn that p + 3p < holds for accelerating expansion. 

Via Hubble parameter, Fredmann equation could be rewritten as: 

k p 



1 (1.22) 



and therefore the cosmic curvature completely depends on the total energy density, 
which is in casuality agreement with general relativity. If the density is above ^Q,k 
will be positive, and results in a limited and closed cosmos; If the density is identical to 



13 



^Q,k will be null, and results in a infinite and flat cosmos; If the density is below ^Q,k 
will be negative, and results in a infinite and open cosmos. So in further discussion it is 
proper to define the unitary quantity which share the dimension of energy density, 
to be the critical energy density of the universe: 

Pcrit = (1-23) 
together with a non-dimensional quantity of relative energy density: 

n = -(^ (1.24) 

Pcrit 

According to the current theories, it is pcrit and ^ that determine the ultimate fate of 
the universe. Hence Friedmann equation take the further form: 

^ Q-1 (1.25) 



thus a clear graph is finished to reflect how the evolution of the universe depends on 
its property parameters, which prepares a firm foundation for cosmic observations and 
cosmological dynamics. 

1.2.4 EOS in Perfect Fluid Scenario 

Actually to resolve cosmological problems under perfect fluid scenario, we still need 
the equation of state (EOS): 

P = Pip) (1-26) 

Thus theoretically, the combination of Friedmann Equation (|1.17p . energy-momenta 
conservation equation ()1.20p and EOS ()1.26p is equal to solve the dynamics of cosmos. 

Generally, a simplest EOS is employed for cosmological perfect fluid: 

p = wp (1.27) 

where G M. For ordinary matter, w = 0. And the accelerating expansion conditions 
p + 3p < can be re-expressed via EOS parameter: 

w<-l/3 (1.28) 

Integration by ()1.12p give rise to the relation of energy density and cosmic scalar factor: 

poca-3(^+'^) (1.29) 
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When the universe is dominated by relativistic matter and radiation (i.e. the baby era), 
w = 1/3, hence 

a(i) cx - ^o)^/^ /9(xa-^ (1.30) 
When the universe is dominated by non-relativistic matter, w = Q: 

a{t) <x. {t - tof/^ , pcxa-^ (1.31) 

All in all, different eras correspond to different domination and evolution. 
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Chapter 2 

Accelerating Expansion of the 
Universe 

2.1 Accelerating Expansion and Dark Energy 

The inflation theory, which is a great success as supplement to Big Bang cosmology, 
declares that the cosmos is flat in early 1980s. Theorists' insistence lead to annoying 
disagreement. However, cosmic observations found that there is far less matter than 
needed for flat topology. Shortly afterwards in early 1990s, COBE is sent by NASA 
to explore the cosmic microwave background. To theorists' delight, the location of the 
first peak in angular power spectrum confirms the flatness of the cosmos. This strongly 
require a new density component to fit the remarkable gap between the matter density 
and the critical density. Before physists could take a calm breath, the redshift mea- 
surement of la supernova declared the accelerating expansion of the universe! This 
milestone is among the greatest ever discoveries, and bring the discoverers the Shaw 
Awards for Astronomy. 

So, CMB exploring proves the flatness, homogeneity and isotropy of the universe, 
and SNlo confirms the accelerating expansion. These results lead to nondebatable con- 
clusions that, if general relativity and the standard model of particle physics are correct, 
then 

(1) The density of cosmos is Q = 1. Since baryons, non-baryonic CDMs and radia- 
tion add up to only around + = 0.3, the major energy (the gap) comprise of an 
unknown component, dubbed Dark Energy (DE). 

(2) Only limited dynamical characters are known. Dark energy has positive energy, 
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negative pressure, and anti-gravitation effects that drive the accelerating expansion. 
DE doesnot condensate at scales smaller than lOOMpc; DE density is almost constant; 
DE is only dominant recently (z G (0 ~ 1) and zq = today), so it doesn't bother the 
nuclear-synthesis process at early era. 

And even up till now, we are still quite unfamiliar with DE. DE may be scalar fields, 
vector fields, tensor fields, or even cosmic strings, etc[7]. And a set of typical data is 
that the critical density is made up of 72% DE, 24% CDM, 4% ordinary matter. 

2.2 Characteristic Parameters of Cosmos 

The properties of the universe need to be parameterized for quantitative study (this 
is a common method to handle huge and complicated systems). 

Presently, there are 12 parameters that are mostly popular. [40]: 

P = {ub, UJc, JC, H,T, Wo, wi,'^mt,, Hs, As, a^, r) 

where tOh = ^l^h? and Wc = O.ch'^ represent the relative energy density of baryonic matter 
and cold dark matter. /C represents the presentday scalar curvature ,and cosmic ob- 
servations support the result of /C = ,a flat cosmosi? represents Hubble parameter;r 
represents luminous depth^ mjy represents the total mass of the three generation of 
neutrinos. . In general computing process, the EOS parameter is often rewritten as 
w = wq + 'Wi{l — a),wo and wi represent the new EOSs respectively. The last foue pa- 
rameters are related to inflation theory.A<j represents the amplitudes of original powern^ 
represents the power index for scalar spectraa<j represents shift of spectra index Ug r is 
the tensor power parameter. 

It is of great importance to investigate these patameters , which help constraint 
the properties of dark energy. Here below we list the data from Willkinson Microwave 
Anisotropy Probe, which can de downloaded from the WMAP homepage of NASA. 

2.3 Mechanism for Accelerating Expansion 

Ever since the proof that DE is dominant, its components, origin, dynamics are 
waiting to be determined. With the parameter constraining via cosmic observations, 
various models have been put forward for a cosmos with accelerating expansion. I would 
divide these models into three categories for further discussion. 
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Table: Data from WMAP. 



Cosmic Parameters 


Symbol 


Value 


-|- uncertainty 


-uncertainty 


total density 


^tot 


1.02 


0.02 


0.02 


Quintessence EOS 


u 


< -0.78 


95%CL 




DE density 


Qa 


0.73 


0.04 


0.04 


Baryon Energy Density 




0.044 


0.004 


0.004 


Baryon Number Dcnsity(cm~'^) 




2.5 X 10-"^ 


0.1 X lO"'^ 


0.1 X lO"'^ 


Matter Energy Density 




0.27 


0.04 


0.04 


Neutrino Density 




< 0.0076 


95%CL 




CMB Temperature(K) 


^cmh 


0.275 


0.002 


0.002 


Baryon-photon Number 


V 


6.1 X 10-^° 


0.3 X 10"^° 


0.2 X 10"^° 


8/i~^MpcSpherical Fluctuation. 


Ss 


0.84 


0.04 


0.04 


Scalar Spectra Indexfco = 0.05Mpc~^) 


Us 


0.93 


0.03 


0.03 


Index shift(A;o = OmMpc'^) 


dug/dlnk 


-0.031 


0.016 


0.018 


Scalar-Tensor Ratio(A;o = 0.002Mpc"i) 


r 


<90% 


95%CL 




decoupling redshift 


^dec 


1089 


1 


1 


decoupling width(FWHM) 




195 


2 


2 


Hubble parameter (reduced) 


h 


0.71 


0.04 


0.03 


Cosmos Age(Gyr) 


to 


13.7 


0.02 


0.02 


Decoupling time(Kyr) 


tdec 


379 


8 


7 


Heavy Ionic time(Myr,95%CL) 


tr 


180 


220 


2580 


decoupling interval(Kyr) 


^tdec 


118 


3 


2 


equality redshift 


^eq 


3233 


194 


210 


Baryonic luminous depth 


T 


0.17 


0.04 


0.04 



2.3.1 Scalar Fields 



Bulk of these proposals are some kinds of scalar fields with shallow potentials and 
thus tiny mass. The field rolls slowly, up (canonical) or down (noncanonical), in its 
own potential. This process takes place only considerable on Hubble time scale. EOS 
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in scalar field models takes the traditionally most simple form p = wp, and EOS pa- 
rameter w turns the identity of different fields. t(;>— lor— l<u;<— 1/3 means 
Quintessence field; i/; = — 1 is Einstein's cosmological constant A; w < —1 reflects 
Phantom field. These cases are all real fields; and those with w crossing —1 can be real 
(Quintom) or complex (Hessence). 

And the key point is that a scalar field is itself meaningful at all, if not treated 
together with the accelerating expansion. Firstly, it's generated or associates with 
many attractive theories, such as superstring quantum field, pseudo-Nambu-Goldstone 
Bosonic Formulation, Brans-Dicke theory. So it's a natural trial to test the cosmological 
consequence of a scalar field, typical work like Ratra and Peebles in late 1980s. Sec- 
ondly, a scalar field could help make clear the CDM spectra. Last but not least, even a 
light scalar field can generate testable consequences in the standard CDM scenario. For 
those reasons, scalar fields such as Quintessence has been discussed pioneeringly, after 
the inflation theory, before the accelerating expansion observations. More discussions 
are available in this thesis afterwards. 

2.3.2 Modified Gravitational Field Equation/Priedmann Equation 

A most remarkable candidate is Cardassian models, which add a new energy term to 
the right side of Friedmann equation, or to the energy-momenta tensor of Gravitational 
Field Equation. This Cardassian term depends on the density of matter and radiation 
only, and the modified Friedmann equation could produce a flat, matter dominated uni- 
verse with accelerating expansion at late time, while the early evolution of baby cosmos 
is not perturbed. The creative work of this thesis is exactly connected with this models. 

Modifications to gravitation equation includes operations on the geometric term of 
Ricci tensor, or energy-momenta tensor. Besides Cardassian, other examples include 
Randall-Sundrum model, etc. 

2.3.3 Modern QFT approaches 

Such models include k-essence, tachyon and Chaplygin gas. K-essence focuses 
on modifying the kinetic energy term of Lagrangian, who acts as the first principle. 
Tachyon plays an important role in standard model of particle physics and superstring 
guage theory, and its connections with accelerating expansion is established after 1998. 
Both k-essence and tachyon modeling are deeply rooted in quantum field theory. 
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Chaplygin gas seems less mathematically complicated than k-essence and tachyon 
formulation. It is a kind of non-relativistic perfect fluid in superstring theory, being 
the only liquid allowing supersymmetric generalization. Its EOS is quite different from 
that of scalar fields, and can be reached via several approaches, such as non-relativistic 
approximation of Born-Inficld theory. In (generalized) Chaplygin gas cosmology, the 
ultimate fate of accelerating expansion is dc Sitter Universe. A most interesting prop- 
erty is that (generalized) Chaplygin gas indicates a unified description of dark energy 
with dark matter. 

In the following test, I firstly review some typical existing work, including four scalar 
fields (Quintessence, Phantom, Quintom, Hessence), k-essence and Chaplygin gas. I put 
my work in the last two chapters, where I come up with a new Cardassian model and 
introduce virial collapse to Cardassian cosmology. 
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Chapter 3 

Quintessence &; Phantom 



3.1 Quintessence Scalar Field 

Quintessence scalar field is the first introduced dark energy model, and has received 
most attention. At early 1980s, a scalar field is employed in the inflation theory, which 
encourages people to analyze the consequence of a scalar field rolling in its own potential 
at the whole Hubble time scales, rather than the tiny time interval at the extremely 
early universe. This work is done in late 1980s, and is re-introduced as the first non-A 
mechanism for dark energy after 1998. 

3.1.1 Basic Structure of Quintessence 

Generally Quintessence is treated as a classical scalar field, with smallest coupling 
with gravitation field, and its action is: 

CcpV-gd^x (3.1) 

where g is the determinant of the matrix g^y. With metric (H ),phenomenently,the 

Lagrangian density with positive energy and positive pressure £^ is: 

^0 = Ig^'d^^d^cb - V{cl)) (3.2) 



Take variance operation on (13. ip . and one gets the energy- momenta tensor for a scalar 
field: 



Tt = d^<i)dy(t) - g^y 



(3.3) 



Take FRW metric (ll.3p in a fiat space into account, via the conservation equation 
T^'^ ~ 0, one can get the dynamical equation for the scalar field </>: 

ij) + 3H<j) + = (3.4) 
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Assuming that the scalar field is homogeneous, whose derivative to spatial coordinates 
is zero and dependent on cosmological time only. Via FRW hne elements in such a field, 
one can get the energy density and pressure of Quintessence field: 

1 :^ 



and 



P4> = i^r + v{^) 



Hence EOS parameter of Quintessence field is: 



(3.5) 
(3.6) 

(3.7) 



The dynamical equations for the evolutionary cosmos is 



a 
a 



~3~ 

SvrG 



02 _ y(<^) 



(3.^ 



(3.9) 



From Eq. (j3.9p indicates that accelerating expansion is possible if cj)^ < V{(j)) holds. Insert 
(13. 7p and ()1.20p into Quintessence field, then: 



= poexp 



/, da 
3{l + w^) — 



(3.10) 



where po is an integral constant. According to ()3.9p . the value of it;^ must be mamong 
-1 to 1. If the rolling in the potential proceeds very slowly, i.e. cp'^ <^ V{(p),th.en 
oj^ = —1 and being a constant; If the rolling in the potential proceeds very sharply, 
i.e. (jP' » y ((/)), then uj^ = 1 and p oc a^^As to the media cases. Quintessence energy 
density will be 

P0 oc a""", < m < 6 (3.11) 

Since u = —1/3 is the critical condition for accelerating expansion and decelerating 
expansion, constraints can be strengthened for < m < 2 for accelerating expansion. 

3.1.2 Quintessence Potential 

In discussions with Quintessence field, three potentials are most popular 

(l)Inverse power style potential by supersymmetry models [lO]: 

V{<t^) = Vo{^Y ao>l (3.12) 
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(2) Exponential potential by Kaluza-Klein theory [TT]: 

V{(t)) = Voe'^-^o (3.13) 

(3) Cosine potential by Pseudo-Nambu-Goldstone Bosonic Mechanism [12]: 

(3.14) 



cos C J ) + 1 



Let's take a closer look through a concrete sample: For power style expansion, the 
driven potential is exponential. The scalar factor expands in power style, a{t) cx t^, 
where p = 1 is the critical condition to separate accelerating and decelerating expansion. 
According to (11.180 . one has 

H = -AttG^'^ (3.15) 
So one could employ the observables H and H to represent the potential V{(p) and 4>: 

This operation can be treated as the reconstruction of the model. Hence the driven 
potential for a power style expanding potential is: 



where Vq is a constant, unitary parameter mpi being Plank mass. In a word, an exponen- 
tial Quintessence potential will drive the universe expanding in power style. Moreover, 
exponential Quintessence could lead to scaling solutions of dynamical phase system, 
where the ratio of Quintessence energy density p,f, to that of background matter pm is 
constant. 



3.1.3 Scalar Fields' problems and Tracker Field 

If one employs a slowly rolling and light scalar field to drive the accelerating expan- 
sion, two serious problem will arise. This is the difficulty that all scalar field come across. 

(l)Fine-tuning Problem 
The kernel of this problem lies that, why the missing energy density differs so much in 
comparison with the typical energy scale in particle physics ? According to thee current 
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data, Q 



0.3, then dark energy density is in the mag nitude of W^'^GeV^, which 



needs an extra 14 magnitude beyond electromagnetic interactions. 

(2) Coincidence Problem 
The kernel of this problem lies at the initial conditions of the cosmos. The evolutionary 
rate of Quintessce field differs very much with that of energy density, and the process is 
quite complicated. So, to evolve into the state where the energy density of Quintessence 
field and matter are of the same magnitude after 137Gyr, the initial conditions must 
be exactly set. 

All scalar fields cosmological models are upset with the two problems, but, the 
discovery of the so-called Tracker Field, which originate from Quintessence and then 
spread to all other scalar fields, throws light on beating these problems|15|. It's found 
that, as for some special potential, such as V ^ Vo</'~") ^ ~ Voexp(M/(/) — 1), F ~ 
Vo[sinh{ay/koA^)]^ [l6j, there are tracker solution for the Quintessence scalar field. 
Tracker solutions have quite unrestricted initial conditions, and the energy density 
ratio of Quintessence to matter is permitted to fluctuate for 100 magnitudes to have 
the same evolutionary results. This means it's almost independent of initial conditions, 
or perfectly Markovian. This solves the two problems referred to above in some sense. 
And from then on, the existence of tracker solution, in the differential autonomous 
dynamical system turn an important judgement for the quality of a scalar field for dark 
energy, and for the quality of a potential for certain scalar field. 

3.1.4 Typical discussions for Quintessence 

(l)In [17] it's shown that the gravitational field equation and conservation laws put 
strict constraints to Quintessence with w = const., For constant EOS parameter w and 
w ^ —1 (simply the case of cosmological constant A): 



where V is the scalar potential for Quintessence field, V' = ±3i7o\/ (1 — '"^^)Vo/2,Vo 
being the presentday value. 

(2)Carroll made typical evaluation to the strength of Quintessence field[18j. Cos- 
mological observations restrict that the scalar potential is quite shallow, which means 
the tiny mass in excited states, = y^V"{(j))/2 < //q ~ 10^'^^ey. In order to support 
the observed energy density, the current value of the potential must be approximately 




(3.19) 
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the closure density, l/((/)o) ~ (W-^eV)^, and the field be (/)o ~ lO^^GeV ~ Mpi = 
(STrG)""*^/^, where Mpi is the reduced Plank mass. A light scalar field will generate 
testable long-range force, and if Quintessence could couple with matter, it would lead 
to the evolution of physical constants at Hubble time scale. 

(3)|19j investigates the origin of Quintessence from breakdown of supersymmetry. 
|20j investigated the possibility that axions and QCD massless quark comprise the 
Quintessence. As a matter of fact, in the following years since the discovery of dark 
energy, attentions are focused on the theoretical formulation, and practical constrains 
for Quintessence takes fall behind obviously. Various connections with Quitessence from 
other branches have been established, and large amounts of theses are available, but 
this thesis wouldnot go further and wider review for Quintessence. 

3.2 Phantom Scalar Field 

According to Riess et al, if the matter density is treated as = 0.27 it 0.04, 
one would see w = — 1.02lQ }g, and —1.46 < w < —0.78 for 95% confidence level. 
If one goes further to take the results of CMB and 2dFGRs into consideration, it's 
w = — l.OS"^'^'''— 0.18. So, w < —1 is also possible. Caldwell for the first time dis- 
cussed the cosmological consequence for w < —1, and introduced the second scalar 
field, Phantom[2T], as opposed to Quintessence. Phantom has many amazing proper- 
ties. Phantom violates the weak energy conditions of gravitation; w < —1 leads to the 
effective sound speed v = \J\d'pldp\ beyond light velocity. The most debating result is 
the so-called Big Rip. The dark energy density is dynamical and increase along with the 
accelerating expansion; thus the increasing dark energy and the accelerating expansion 
unite to form a positive feedback pair, and one day in the future, the huge negative 
pressure will tear out the spacetime and all matter. For a review, see [22] and [2l] . 

3.2.1 basic Structures of Phantom 

Canonical scalar fields couldnot realize the process for w < — 1, so noncanonical ones 
are employed for Phantom. A noncanonical Lagrangian for free Phantom field, which 
also has the smallest couplings with gravitation, can be written as|21][25]: 

ll^ = -\9^''d^,4>du4>-V{4>) (3.20) 

Hence the dynamical equation for Phantom field is: 



+ mi\) - 




(3.21) 
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the minus sign before thee potential term makes Phantom climbs up the potential, as 
opposed to a canonical field's rolling down in its potential (see Quintessence). 



Variance of Phantom Lagrangian leads to Phantom energy density and pressure: 

(3.22) 



and 



P4> = -\^''-yi<i>) 



(3.23) 



as well as Phantom EOS (To avoid the increase of Tackyon module, w is always assumed 
to be constant in Phantom): 



W4> 



P4> 



(3.24) 



Obviously W(f, < —1 when 5^ < ^(0)- -^^d energy-momenta tensor for Phantom field 
is: 



^/t = -dn<pd^(f> - guv 



(3.25) 



3.2.2 Hubble Parameter &; Phantom Cosmology 



The energy density of matter becomes equal with that of Phantom field at the 
cosmological time teq-, and the scalar factor driven by Phantom energy will be: 



a{t) ~ a{t 



eqj 



(1 + w) w 



3(1+10) 



w <-l 



(3.26) 



So the scalar factor will be divergent one day, i.e.t — ^ tsR = ^1^^ ieg and a{t) 
Insert w < —1 into Friedmann Equation: 



00. 



H{z) = Hq + zf + ^x{l + z 



,3(l+w) 



1/2 



(3.27) 



where Hq = H{z = 0) means the presentday value of Hubble parameter, and w = wde 
refers to EOS parameter of dark energy. Q„ 



^iGporu^ ^ s^GpuE ^ is relative energy 



density for matter and dark energy respectively, and Jl^ + = 1- Thus, 

-1 



H{z) - t-^ 



t 



(1 + w) w 



(3.28) 



When t tsR = ( teq-, the Hubble parameter is also divergent mathematically. 



This indicates that one day the expanding velocity will approach the extremely upper 
limit. So Phantom dark energy density p{t) oc [$7^(1 + 2;)^ + 0^(1 + 2)^*^^"'""')] ^ is 
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divergent at Big Rip, too. At the Big Rip singularity, spacetime and matter are all 
destroyed. 

Phantom is also a kind of scalar field under large amounts of discussions, and its 
connections with quantum field theory, standard model of particle physics, superstring 
theory and quantum cosmology have all been established [H3] [HI] [35] [HH], etc. And 
not all Phantom cosmos of ii;^ < — 1 ends up with Big Rip, but extra mechanism must 
be introduced, see [26] [27] [28] [29] [30] [H] [32] and the pioneering work of classification 
of Big Rips from Qhantom and Quintessence theories [23], as well as alternative attempt 
to quintom [27] . 
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Chapter 4 

Quintom Hessence 



4.1 Quintom Scalar Field 

4.1.1 Scalar Field with w crossing -1 

Riess et al analysed their SNla Gold data (with dark energy EOS parameter w{z) = 
wo + w'z), and found w = — 1.021q^9, which indicates that the presentday EOS param- 
eter is probably crossing -1. Shortly, Alam, Sahni ct al combined Gold data with CMB 
data to examine the behaviors of dark energy, and they drew clear conclusion that EOS 
parameter evolves continuously, from iy~0at2;~l, tow^— latzciO today. 

Hubble parameter and luminal distance is connected via: 

Alam et al selected a formalism in which Hubble parameter is independent of models: 

H{x) = Ho[nmox^ + Ao + Aix + Asx^]^ (4.2) 

where x = 1 + z, and ^Imo being the presentday matter density, and Aq, Ai, A2 are 
constants. Thus EOS parameter is reconstructed tobe: 

(2x/3)H'/H -1 

where prime ' refers to derivative with x. Together with modified test, the above 
conclusions are drawn. 



Hunterer confirmed the result of crossing -1, via data fitting of piecewise function 
and cosmic obscvations, and to be more exact, the cosmological redshift z < 0.2 for 
w < —1, and z > 0.2 for z > —1. Corasaniti went even further to make combined 
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constraints with CMB, SNla, and large scale exploring data, to get the result of w 
crossing -1. 



As it is put above, —l<w<—^ for Quintessence and w < —1 for Phantom, and 
this constraints is principly inviolated, so, both Quintessence and Phantom are power- 
less to realize the crossing of the state w = —1. And therefore, the so-called Quintom 
scalar field is put forward, which aims at crossing EOS -1. 

However, to build such a Quintom model, one runs into the No- Go theorem (see A. 
Vikman|37j): (l)Under the frame of four dimensional FRW metric and general relativ- 
ity, (2)if dark energy is a single perfect fluid, or single scalar field (p with Lagrangian 
C = C{(j), d^(j)d^(j)), (3)with smallest coupling to gravitation, then the EOS parameter 
w couldnot cross -1. This means, to make w cross -1, at least one of the three condi- 
tions should be violated, such as Quintom with double field [38j. or Quintom based on 
Lagrangian with high order modification term [39]. 

4.1.2 Basic Structure of Double Field Quintom 

Quintom field includes two real scalar fields, one being canonical Quintessence style 
field (j) )the other being Phantom style field 93. 

Lagrangian deensity of Pantom field is: 

C = ^d^^Qd^4>Q - ^d''(t>pdf,(Pp - V{(t>Q, (pp) (4.4) 
thus Quitom dark energy density is: 

] 2 ■ 2 

/'9 = y- Y + mv') (4.5) 

and with pressure of: 

] 2 • 2 

pq = \-\-y{4>,^) (4.6) 

With the dark energy density and pressure, EOS parameter of Phantom dark energy 
is: 

Pq <P^ -ip^ + 2V{(P,ip) 

From Eq ()4.7l) we can see if <j)'^ > ip'^, u; > —1; if c/)^ < (/9^,cj < —1. EOS of Quintom 
crosses -1. 
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At first these two componential potential have no interactions with each other, 
Quintom potential is the simple addition of those two potentials: 



V{<t),^) = Vi{(t>) + V2{^) 



(4.8) 



take a simple potential V{4)Q^(t)p) for example: 

<t>p) = exp{-^4)Q) + exp{ 



Mp 



A 



(4.9) 



For these cases, the steady attractor solution is in Phantom style, or the ultimate state 
is Phantom, where one has to face the nightmare of Big Rip. 

4.1.3 Generalized Quintom 

There are various methods to make expansion, a simplest and natural way is to set 
weights for the addition of Quintessence and Phantom. 

4.2 Hessence Scalar Field 

Quintessence, Phantom, Quintom are all real, and the earliest attempt to describe 
dark energy is Spintessence model. However, Spintessence falls into deadly trouble 
when dealing with the steadiness of Q-balls and is generally abandoned as a dark energy 
model(but Spintessence turns out to be a good candidate for dark matter). 

4.2.1 Basic Structure of Hessence 

Considering a single complex scalar field with EOS parameter crossing -l|41j[32]: 



But, Dark energy based on canonical Lagrangian is not steady (Spintessence). Thus, 
the modified noncanonical Lagrangian is introduced as the first principle for Hessence 
complex model: 



<I> = 01 + i4>2 



(4.10) 



As for complex field, its canonical Lagrangian is: 







(4.12) 



the action of Hessence field is: 




(4.13) 
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where {tp, 0) are two newly defined variables for Hessence: 

01 = (/)cosh6', (/)2 = </) sinh 6* (4-14) 

thus 

cp'^ = - (Pi coth0 = ^ (4.15) 

02 

This is why Hessence receives its name, and H here stands for hyperbolic. 



In flat FRW spacetime with scalar factor a{t), and supposing (p and 9 is isotropic, 
from Hessence Lagrangian (j4.12p . one gets Hessence dynamical equation which depends 
on (j) and 9: 

dV 

(P + 3H(P + (P9^ + — = (4.16) 
o(p 

^'^9 + {2^(P + 3H(P^)9 = (4.17) 
Thus Hessence energy density and pressure are: 

Phe = ^(<A'-</''^') + n^) (4.18) 

Phe = l{^^-(t>'e^)-V{^) (4.19) 

Based on (l4.17|) .there's 

Q = a^\i!^9 = const. (4.20) 

where Q is the conservational quantity in certain volume, which is generated due to 
inner symmetry of complex scalar fields. In details, the Hessence Lagrangian (|4.12|) 
is conserving under local gauge transformations, or remain a same mathematical form 
after (p ^ (p and 9^9 — ia. According to Noether theorem, conserving fluid and quan- 
tity associates (noted Q). On the other hand, Hessence potential V{^, <!>*) or V{(pi, (p2) 
should only depend on + or cpf — (j)2{ov noted V{cl))), to keep the symmetry 
holding. Now turn attention to the comparison of Hessence Lagrangian with Quintom 
Lagrangian, which shows that it is 9 term or Q term that plays the role of Phantom 
field. When Q = 0, Hessence reduces to Quintessence, and there is no more necessity 
for crossing w = —1. And there's the pity that physical meaning of Q is still uncertain. 



It should be stressed that Hessence differs from Quintom and is an independent 
model. An important difference is that (pf — (p\ is treated as a whole part to be the 
independent variable, while the two componential fields are generally independent to 
each other in Phantom. 
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With Eq ()4.20p one can rewrite the dynamical equation, energy density and pressure 
of Hessence as: 

dV 

I t is notable that (|4.2ip is equal to Hessence energy conserving equation: 

Phe + 3H{phe+Phe) = (4.24) 
while Friedmann and Raychaudhuri equations turn to be 

H' = ^{phe + Pm) (4.25) 
H = -AirG^Phe + Pm+ Phe + Pm) (4.26) 

where pm and pm are the energy density and pressure for background matter. Ulti- 
mately, EOS parameter of Hessence is: 

And apparently when (p^ ^ l^W^i Hessence EOS parameter Whe ^ — Iwhen (p^ ^ '^W^i 
Whe ^ — 1- Via the evolution of cj) field, EOS parameter w crosses -1 when c^^ ^ 
Eqs I4.20M.271 all depends on straightly rather than Q, which means +Q and 
—Q share the same dynamical behaviors, and leads to the assumption of existence 
of Hessence DE and anti-Hessence DE. 



Canonical and complex fields could easily generate the so-called Q — hall structure. 
Q-ball is a kind of non-topological soliton, whose steadiness is related to the corre- 
sponding Q conserving quantity. As to complex fields with canonical Lagrangian, such 
as Spintessence, once the energy density fluctuates, the perturbation would magnify 
nonlinearly at an exponential speed to form Q-balls condensation. This a a big trouble 
for DE, because DE doesn't condensate at the scales smaller than lOOMpc. Once come 
into being, Q-balls act like matter and its energy density decay into other particles 
at the speed proportional to a~^. However, as for Hessence DE, whose Lagrangian 
is noncanonical, the formation of Q-balls are easily avoided. Although the Hessence 
Lagrangian never occures before in quantum field theory, due to its steadiness conse- 
quences, Hessence is treated highly as a reasonable complex scalar field candidate for 
dark energy [35] [36] . 
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4.2.2 Hessence Reconstruction via Hubble Parameter 

Firstly,denote M^^ = {SnGy^^^, and M^^ is the reduced Plank mass. Then ([32S1), 
(j4.26p could be rewritten as [33]: 

= ^^(Phe + Pm) (4.28) 

H = -^{Phe + Pm+Phe+Pm) (4.29) 

This is an interesting formulation. According to the two equations above, one gets: 

1 

2 



V{^) = SM^iH^ + M^iH - (4.30) 



^ = -2MliH-pm (4.31) 



and the redsliift being z = a ^ — 1 (supposing zq = 1 and subscript refers to today), 
then for arbitrary /, 

f = -{l + z)Hf (4.32) 
az 

and (j4.3Up . (|4.3ip could be expressed as: 

AH 1 
"dJ ~ 2 



V{z) = SM^iH^ - M^iil + z)H^ - ^/5^o(l + zf (4.33) 



) ' - ^(1 + z^H-' = (1 + z)-^hJ-^ - pMI + z)H-^ (4.34) 



now introduce such non-dimensional quantities as below: 



^=1^, Q = QMliHo (4.36) 
-no 



thus (051) . (fOil) turn to: 



dfJ 3 

y(z) = 3^2 - (1 + z)i?— - -J7^o(l + ^)' (4.37) 

f ? I - Q''^"'(l + ^)^^"' = 2(1 + z)-^H^ - 3J7^o(l + z)H-^ (4.38) 
\dz I az 

where QmO = Po/i^^piH^) refers to the presentday energy density of matter. Once 
H{z) or H{z) is known, we could reconstruct V{z) and </)(z) via (j4.37p and (j4.38p . and 
reconstruct V{(l}) via V{z) and (/'(z). From (I4.3ip (j4.37p one gets the reconstruction of 
Hessence EOS parameter: 

wUz) - ^ = "' + ^^1 + "^^ (4.39) 
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Kinetic Energy could be induced from reconstruction of ()4.3ip : 



a2 ml 



K = ^- (4.40) 



thus, 

And this reconstruction depends on testable H{z) rather than certain models 



K{z) = = (1 + ^)iH)^ - ^n^oil + zf (4.41) 



4.2.3 Autonomous Equations for Hessence 

Firstly, to express the background matter via EOS of perfect fluid: 

Pm = WmPm = {w - l)pm (4.42) 

where < 7 < 2. Particularly, 7 = 1 and 7 = 4/3 represent the matter dominated 
state and radiation dominated state respectively. Introduce the interaction term C of 
Hessence and matter at energy equilibrium, and when energy conservation conditon 
ptot + 3H{ptot +ptot) = holds 02] [m, there's: 

Phe + 3H{phe + Phe) = -C (4.43) 
p„, + m{pm +Pm)=C (4.44) 

If there's no interaction between Hessence and background matter, C = 0, and interac- 
tion for C / 0. 

To establish the differential autonomous equations, one takes the following trans- 
formations, = 87rG,K > 0, and introduce the non-dimensional variables: 

Kcj) kVV k^/p^ ^/6 K Q 

X = —;= — , y = —pz — , z = — 1= — , u = — -, V = —= 5"- (4.45) 

y/6H y/3H V3H nek VQH a3<l) 

They forms the Hessence phase space, and now Hessence dynamical problems can be 
treated as geometrical problems in Hessence phase space. And according to basic 
Hessence equations one gets 

„2 „,2 , 7^2 T\ „.„.2 /^„.2. 



x' = 3x{x^ -v^ + -^z^ -l)-uv^ -J -y^f - Ci (4.46) 



y' = 3y(x2 -v^ + Iz^) + ^ ^xyf (4.47) 
z' = 3z{x^ -v'^ + Jz^ _ 2) + (272 (4.48) 
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u' = -xu^ (4.49) 
v' = 2,v{x^ -v^ + - 1) - xuv (4.50) 
where prime operator (') refers to derivative to e-folding time N = Ina, and 

Ci ^ C2 ^ = -Ci (4.52) 

Thus Friedmann Equation reads: 

+ y'^ + -v"^ = 1 (4.53) 
And the relative energy density for each component reads: 

Q,, = ^=x2 + y2_^2^ ^^=PlI^ = ,^ (4.54) 

Pcrit Pcrit 

and the effective EOS parameters for Hessence and the whole cosmos are: 

Phe x^ -v'^ - y2 

Whe = = -o 2 y^-^^l 

Phe x^ -v^ + 

Phe~\~Pm 2 2 2,/ i\2 r-c\ 

Weff = = x —V —y +(7 — 1)2; (4.56) 

Phe + Pm 

the critical points(x, y, n, v) is gained via setting x' = y' = z' = u' = v = {). And this 
is the foundation for further analysis. Hessence has perfect mathematical structure, so 
it is employed to show the normal process of dynamical analysis here. 



It is shown in [l2] that, as for exponential potential and power potential, all late 
attractor solutions meet uj^e ^ —1 and uJeff ^ ~1 with 4- form interaction term C. 
Nevertheless, Phantom style solutions with Lv^e < — 1 or cje// < — 1 is not stable, 
and will evolve into Quintessence attractor with ujhe > — 1 or ooeff ^ —1; or de Sit- 
ter attractor with uj^^. = — 1 or / = ~ 1 • So Big Rip is avoided in Hessence cosmology. 

In [U], it is shown that in 4 dimensional phase space, the existence of stable late 
attractor solutions is independent of the forms of ¥{(()) and C, but the stability of the 
solutions depends on the second order derivative of Hessence potential V{4)). 
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Chapter 5 



k-essence &; Chaplygin gas 



5.1 k-essence Mechanism 

The Scalar field dark energy models employ slowly rolling potentials to drive the 
accelerating expansion. Yet, such a cosmos scenario can be gained via modifying the 
kinetic energy term. In 1999, Armendariz-Picon realized that the kinetic energy term 
could lead to cosmic inflation in the baby era that is governed by high energy mech- 
anism, and introduced k-inflation. Afterwards Chiba utilized the same mechanism to 
describe dark energy that drive the accelerating expansion, and Armendariz-Picon im- 
proved this work into k-essence. Fork-essence, high order modifications appears at the 
kinetic terms to drive the acceleration. And Bohn-Infeld field is one form of k-essence. 
If the form of pressure is properly selected, the fine-tuning problem and coincidence 
problem can also be solved [17] [H] [l9] . 

k-essence has noncanonical kinetic energy. As for a scalar field with (j) and kinetic 
energy X = — (l/2)(Vi;^)^, the action with the most generalized form is: 



where Lagrangian density is the same with pressure density p{(j),X). and the action is 
also used in generalized Quintessence models. 

Usually, Lagrangian density of k-essence obeys the following constraints: 



This constraints originates from superstring theory, in which low energy effect will leads 
to higher order derivative terms from a' and loop modifications (the relation of a' and 





p{4>,x) = f{4>)p{x) 




36 



string length scalar As is a' = \s/2tt). The four dimensional effective string action is: 
S = j d''x^g{Bg{ct>)R + BfmV<pf-a\c'~l^B^l\ct^){V<^^^ (5.3) 

where A;^ = SttG = 1. Here cf) is Dilaton field, which controls the string coupling 
strength g1 via gl = . For weak couplings « 1, the coupling equation is Eg 
B^^'^ ~ B^^'^ ~ e~'^ . If the coupling is one order, it will take a more complicated form. 
Via conformal transformation g^j^y = Bg{(j))g^y, the string action will transforms into 
Einsteinian action form: 

Se = j (fx^[]^R + K{4>)X + L{4>)X^ + ...] (5.4) 

with the component 

L(0) = 2cf)a'4^)(0) (5.6) 
Thus we could get the Lagrangian with noncanonical kinetic energy term: 

p{<l), X) = K{<^)X + L(0)X2 (5.7) 

and redefine the field (l)new as: 



L 

\K\ 
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Thus Lagrangian transforms into: 

p{^,X)=f{<p){-X + X^) (5.9) 

where (j) = (t)new,X = X^ew = (L/|K|)Xoid,and /((/>) = K'^{4)oid)l^{.4>oid)- This is the 
k-essence model induced from (j5.2p via ]5(X) = —X + X^. Thus the pressure is also 
known, thus the energy density of the scalar field ^ is: 

P = -p = fm-X + 3X2) (5^^Q) 

With energy density and pressure, one immediately gets EOS parameter for (p: 

p 1-X , , 

and as we can see, for certain X, W(f, is also a constant. For example, when X = 1/2 we 
could get a cosmological constant style EOS parameter = —1; to get an accelerating 
expansion w < —1/3, one gets X < 2/3. 
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According to energy density p, one gets the continuity equation ()1.20p . When matter 
and radiation is dominant and EOS parameter for background fluid is Wm, Eq (jl.l8p 
gives rise to the Hubble parameter: 

2 

H = — (5.12) 



and the energy density being 

. _ 2{l + w^) 



(1 + Wm){t - to) 

For certain X, or certain w^f,, f{(j)) is constraint for 



(5.13) 



/(</.)« (0-00)-", a = ^^^^^ (5.14) 

when W(f, = Wm, /(</>) esc (0 — ^o)"^ for radiation dominated or matter dominated 
era, and this agrees with the scaling solutions. Thus, as for dark energy, one must take 
fine tuning of /(</>) to satisfy the presentday cosmic density. And one should note that 
the density of the scalar field must be significantly less than that of the background 
fluid, i.e. p <^ pm, which contradicts with a DE dominated universe. For example, 
f{4>) (X {(f) — 0o)^! except for a late scaling solution, we has another solution of acceler- 
ating expansion. Actually this is the divide of accelerating expansion with decelerating 
expansion. 

Prom (jS.lip we could know that, the kinetic term x plays the key role to influent the 
EOS. If 1/2 < X < 2/3 holds, the scalar field (f) will acts as dark energy with ^ a ^ 2. 
Armendariz-Picon et al have made more generalized analysis with p{X) to avoid the 
fine-tuning problem. 



5.2 Chaplygin gas model 

In the discussions with Quintessence, the background perfect fluid always has a 
EOS of p = wp, and change different potentials and Lagrangian. What if employing 
some other form of EOS ? It is very coincident that, when three superstring experts 
were studying the stability of black holes in brane world, they found they have to add 
a background matter called Chaplygin gas, to ensure the stability. Chaplygin gas is 
a perfect fluid with negative pressure in superstring theory. They proceeded to calcu- 
late the cosmological consequence for a PRW universe comprised of matter, CDM and 
chaplygin gas, and found that a flat universe with accelerating expansion is gained, and 
chaplygin gas not only acts as dark energy, but strongly indicates a unified description 
of dark energy and dark matter. Hence, Chaplygin gas is introduced as a DE candidate 
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in a quite natural manner. 



5.2.1 Special Chaplygin gas 

EOS of Chaplygin gas is (all components depend on comoving coordinates): 

PCG = (5.15) 

PCG 

where pcG > 0, pcG > 0, ^ is a positive constant. Chaplygin gas obeys the energy- 
momenta equation (|1.20p . ()1.12p . Insert the EOS of Chaplygin gas into ()1.12p . and the 
integral induces: 

PCG = + ^ (5-16) 

Where B is an integral constant, and is selected to be positive in the following discus- 
sions; a is the scalar factor. Eq (l5.16p show the behaviors: 



a^{^)^^pGC 3- (5.17) 

a>(^)^ — 'PCG^-PCG^^ (5.18) 

And we can easily draw the conclusions that, at early epoches, a is small enough to 
ensure a <C (§)6 , and density of Chaplygin gas pcG ~ which is like non-relativistic 
matters; lately a is large enough to ensure a 3> ^)6,thus pcG ~ —PCG ~ and acts 
like the cosmological constant, which has negative pressure to drive the accelerating 
expansion, and the time needed to arrive at this state is: 



1 / (M + 4 + ^ / B 

In ^ 2 arctan (/ 1 + -— 1 (5.19) 



Because Chaplygin gas acts like dark matter and dark energy at different epoches, it 
indicates a unification of CDM and DE. 

A more detailed perturbative form of (|5.18p is: 

pr^VA+J^a^ (5.20) 
V 

p~ (5.21) 

(I5.20p . (I5.2ip show that the expanding universe is made up of two components, the 
major part acts like the cosmological constant with p = —p; the minor part acts as 
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p = p, whose energy density decrease sharply with the expansion. 



The cosmological consequences of Chaplygin gas could be treated as a homogeneous 
scalar field (/>(t) togther with a potential V{(j)). Consider the Lagrangian: 

J^GC{4>) = - V{<P) (5.22) 
and set the density of the field equal to that of Chaplygin gas: 



P^ = \^^ + V{cP) = ^Ja + ^ (5.23) 



the pressure is identical to Lagrangian: 



After some algebra, the reasonable potential is: 

VW = ivl(™h30 + ^) (5.25) 

The potential is independent of B, and it reflects the properties EOS of Chaplygin gas 
(f5?T5]) only. 

There are several approaches to EOS of Chaplygin gas, including: 
(l)Action and Potential of Quintessence 

''W = f("="* + ^^' (5.26) 

and in return, with the factors A, B and scale factor a in Chaplygin gas, one can 
re-express the kinetic energy and potential in Quintessence: 

B 



12 



and 



(2)Born-Infeld Lagrangian 



'{A + B/a^ 
2a^y/A + B/a^ 



(5.27) 



y(0) = + (5.28) 



£ = -Voy/l-d^dt^ (5.29) 

has the induction that: 



Vq 



(5.30) 



p = -Vo^/l - df,d>^ (5.31) 
which end up with Chaplygin gas EOS. 
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5.2.2 Generalized Chaplygin gas 

Firstly to generalize EOS of Chaplygin gas(GCG) into: 

p=-^ (5.32) 

where < /3 < 1. (further consraints sre mage in [24] to get ^ /? ^ 0.5. The following 
task is to find out the scalar field corresponding to (|5.32p 

For a complex field with non-zero mass, supposing: 

<I> = ^-^m^ exp {—imO) (5.33) 

with Lagrangian of: 

Cgcg = 9^" {^^0,^9,,) - y(0V2) (5.34) 

where , fi refers to variance with spatial coordinates . The condition (j)^^ <C rrKp is 
used to ensure the fluctuations of energy density in spacetime, and together with the 
generalized EOS (j5.32p and Thomas-Fermi, to rewrite its Lagrangian: 

J^TF = y5^"^,m^,. (5.35) 
This acts as first principle for generalized Chaplygin gas cosmology. 

From (l5.32p . one has a generalized ()5.16p : 

/ B \^ 

Energy density for generalized Chaplygin gas being: 

PGCG = iA + ^^)^ (5.37) 



the effective EOS is: 



a3(i+/3) 

and the perturbations ()5.20p . ()5.2ip turn to be: 

f I \ B 



co{a) = '"^"^ I (5.38) 



1 + 



p ^ + ( — ^ ) ^l^a-=^(^+^) (5.39) 



^ f I \ B 



p ^ _ATTJ^ + ( I ^a-3(i+« (5.40) 

Thus, as for the generalized Chaplygin gas, the expanding cosmos is also made up of 

1 

two components. The major part is dark energy with p = p and . The minor part 
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is of EOS p = (5p now, and the potential for generalized Chaplygin gas with Lagrangian 
(I535D is: 

where ^ = B^^-PI^+P^a^^^-^^^ . 



There's some other generalized form, for example 

Aia) 

p = — 5.42 

P 

where A{a) depends on the scale factor a, and is set positive. A modified or generalized 
version can be more flexible and can supply better cosmic scenario. 

5.2.3 Generalized Chaplygin gas &; Hessence 

If the field is properly selected, the late behaviors of Hessence complex scalar field 
could be expressed via generalized Chaplygin gas. Firstly, Hessence should obey EOS 
([532]): 



Phe = — J- (5.43) 
Phe 



Via EOS of Hessence g^Z]): 



2V{(I)) = Phe - Phe = Phe + ^ (5.44) 

Phe 

6^2 = Pf^<^ + Pf^^ = Phe J- (5.45) 



" f^he 



Friedmann Equation p.l7p for Hessence dominated 



If the above condition holds, 



a 

Analyzable solutions are achieved: 



- =H^ = ^Phe 5.46 
a 6 



GttGBQ 



2 



a^2{A-B/a^) 



(5.49) 
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5.2.4 Chaplygin gas & Superstring QFT 

The origin of Chaplygin gas cosmology has shown its connections with superstring 
quantum field theory. Chaplygin gas is the non-relativistic limit of Born-Infeld mecha- 
nism. For excellent details, see Jackiw[52]. 
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Chapter 6 



Hyperbolic Cosecant Cardassian 
Cosmology 

6.1 Introduction 

In 2002, K. Freese and M.Lewis from Michigan University came up with a new 
mechanism for dark energy, where dark energy doesn't exist at allj53j! They modified 
Friedmann equation, or the gravitational field equation, which is of fundamental sig- 
nificance in cosmological dynamics, to establish a scenario where the universe is flat, 
matter dominated, and has accelerating expansion recently. They treated matter and 
radiation together (denoted as a matter term), and introduced a new term, Cardas- 
sian energy density, which depends quantitatively only on the density of the matter 
term, to the right side of Friedmann equation (the original discussion takes power func- 
tion style Cardassian term as an example). When dated back to the baby era, the 
Friedmann equation would return to its normal form with no appearance of Cardassian 
term, and the process such as nuclearsynthesis and galaxies formation are not disturbed. 
Cardassian term turns dominant only recently, with redshift z ~ 0(1)- In Cardassian 
cosmos, the critical density is remarkably less than the normal value, and can be exactly 
^crit = via tuning parameters. So, in Cardassian world, matter itself determines 
the flatness of the universe, and Cardassian term takes over to drive an accelerating 
expansion recently. There is no room for dark energy any more. Cardassian cosmology 
has survived several observational test, such as the age of the universe, CMB, etc. Ini- 
tially, the origin of Cardassian term is attributed to the natural consequence when we 
insert the observational cosmos into a world with higher dimension. 

But the work of K. Freese and M. Lewis didn't raised much attention originally. 
So, 7 months later Freese expressed the idea again and generalize the density term in 
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Friedmann equation to a general form g{p) by the way [5l]. Shortly, Gondolo and Freese 
put forward the relativistic fluid interpretation of Cardassian dynamics [55] (denoted as 
GF fluid hereafter). Cardassian cosmology is raised to the hight of theory from phe- 
nomenology, and the dynamical mechanism is analysed comprehensively. This work is 
a milestone in the study of Cardassian cosmology, which lay the foundation for future 
research on this subject, and up to now, nearly all the discussions with Cardassian are 
based on GF fluid scenario. To ensure the positiveness of sound speed, they analysed 
a new kind of Cardassian term, the polytropic and a better modified polytropic style. 
They also come up with a proposal that the negative pressure arise from the interactions 
between dark matter. From this thesis on, Cardassian cosmology wins large attention 
and hot discussions. 

In [SB] it is shown that power style Cardassian violates the weak energy conditions 
of gravitation on small scales. Although suffering this problem, power style Cardas- 
sian, due to its simplest form, is still widely employed when new ideas are introduced. 
Two other important work attribute to [57j(M Szydlowski and W Czaja) and [58] (R 
Lazkoz, G Leon). The former provides an excellent example for autonomous dynamics 
and numerical analysis for Cardassian. The latter divides the Cardassian discussion 
of autonomous dynamical system into high energy limit and low energy limit two cat- 
egories for the first time. The typical work to generalize Cardassian model includes 
the introduction of exponential Cardassian term[6T]; Statefinder diagnose for modified 
polytropic Cardassian [62]. Generally, computations with Cardassian cosmology require 
the background to be homogeneous and isotropic, and thesis [63] makes an attempt to 
cancell the restriction of homogeneity. Meantime, numerical analysis and observational 
constraints for Cardassian take step [Mj[65][66][67][M][69][70][7l][72][73][75 
|75j [76j [77], mainly with CMB, the cosmic age, gravitational lensing and large scale 
exploring. 

Freese and her cooperators have written several reviews on Cardassian cosmology, 
which does good for a glance of Cardassian, see [5l] [59j |60| . And in this chapter, we 
will do a job similar to [61j , to introduce a new model via the introduction of hyperbolic 
cosecant Cardassian term. 

6.2 Fundamentals of Cardassian Cosmology 

As expressed above, the following analysis is based on GF fluid scenario. In the orig- 
inal work, matter and radiation are treated together for a matter energy density term 
and the independent variable of Cardassian term. The following analysis will change 
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into two cases. This treatment is reasonable, because Cardassian term is effective only 
recently {z ~ 0(1)) when matter is far more enriching than radiation, and radiation is 
considerable and makes a difference only in the early era, when the Cardassian modi- 
fication is neglected. But for a clearer scenario or as a generalization, to treat matter 
and radiation separately is also a good choice. So in this text two cases are discussed, 
namely treating matter and radiation together as a matter term and treating them 
separately. 

6.2.1 Fundamentals for Implicit Cardassian term 

In Cardassian cosmology, the total effective energy density is|53j [53]|55j : 

p' = g{p) (6.1) 

where p = Pm + P-y reflects the summary of energy density of matter pm and energy 
density of radiation , which are testable parameters. From the first law and second 
law of thermaldynamics, we have Td{sV) = d{pV) + pdV; for adiabatic expanding 
d{sV) = 0; conservation of particles' number d{pmV) = 0; conservation of radiation 

4 

d{pjVs') = 0. Insert these relations into ()6.ip and we could get the total pressure: 

p=-p+p^- — + (6.2 
dpm 3 ' dp^ 

In Cardassian cosmos, Friedmann equation reads: 

= '-^9{P) (6.3) 

and Raychaudhuri equation reads: 

dp' 4 dp' 

the energy-momentum conservation equation: 

p' = -3Hip' + p') (6.5) 
The above equations are the fundamental equations for Cardassian cosmology. 

6.2.2 Fundamentals for Explicit Cardassian term 

Firstly dividing p' = g{p) into two terms: 

p' = P + f{p) (6.6) 
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The first term p is the matter and radiation term, and f{p) = g{p) — p is the exphcit 
Cardassian term, which is the modification term of Friedmann equation. When matter 
and radiation are treated separately: 

p' = Pm + Pi + fiPm, P-y) (6.7) 

Thus from (j6.2p one gets the total energy: 

,' = El + ,^lL + *^^f (6.8) 

3 dpm 3 'dp., 

and from this equation, we can easily see that only radiation and Cardassian term 
contribute to the total pressure, or EOS parameter for matter is if = 0. This conclusion 
will be quoted in next chapter for virilization. Friedmann equation and Raychaudhuri 
equation are 

= —^[Pm + P-y + f(.Pm,P^)] (6.9) 

and 

4 df A df 

i/ = -4vrG(p. + -,, + p„^ + -,,—) (6.10) 

And energy conservation is equal to the conservation of particles' number and conser- 
vation of radiation: 

Pm = -SHpm (6.11) 
Py = -AHp^ (6.12) 

6.2.3 Phase Space, Autonomous Equations for Separate Matter &z 
Radiation 

Based on Friedmann equation, we define 

SttG 

^ = ^Jp[Pm + P'y + f{pm,Py)] (6.13) 

where Clm, ^7, ^card are the relative density parameter of matter, radiation, Cardassian 
term respectively. Now we introduce the non-dimensional parameters (with = SttG): 

where N is e-folding time. Apparently N keeps the same pace with the cosmological 
time in our flat, expanding FRW universe. Those newly introduced variables form the 
4 dimensional phase space, and when x,y,z, N are inserted into the basic dynamical 
equations, they will turn to the geometric curves in the phase space by: 

— = + -x-" + 2x?/ h2xy — (6.15a) 

dN 2 2 2 opm op^ 
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These curves are called the differential autonomous equations of the system. For further 
discussion, we treat the cosmos as flat (the confirmed result by CMB etc), then there 
will be the restriction + = 1, which simplify the equations (j6.15p above into: 

dx 3 3 3 2 3 3 5/ 2 "9/ fata \ 

— = + -x-' + 2xy +-2;^ ^2xy — (6.16a) 

dN 2 2 2 opm op^ 

Particularly, if the Cardassian term has the symmetry fipm-.P-f) = f{Pm + P'y), tbe 
equations (I6.16P can be further reduced into: 

^^4+(2.' + 2./)(i + |L) (o,,) 

+ + (0,7b, 

This case is important, because the hyperbolic cosecant Cardassian term, as we will 
put forward later, just has such a lovely character. 



So, once f{pm,P-y) is known, we will master the dynamics of the Cardassian universe, 
and the proximation behaviors of the corresponding autonomous equations in phase 
space. 

6.2.4 Unitary Treatment of Matter &z Radiation 

If matter and radiation are treated together as a matter term, or only take matter 
into account for recent universe, the energy total density will be 

p' = Pm + f{Pm) (6.18) 

From (|6.2p one gets: 

p' = Pm^-f (6.19) 
OPm 

Friedmann equation and Raychaudhuri are separately: 

H^=^-^[Pm + f{Pm)] (6.20) 

and 

df 

H = -47rGp™(l + -^) (6.21) 

dPm 
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Now the energy conservation equation p' = —3H{p' + p') is equal to the conservation 
of particles' number pm = —3Hpm- When matter and radiation is treated together, 
there's convenient form for the speed of sound, which should be positive: 



2 _9p' _ W' 



ci = ^= Pm^ (6.22) 

dpm 



Based on Friedmann equation, we define 



l=^^[Pm + f{Pn.)] (6.23) 
= ~l~ ^card 

and introduce the new variables {k^ = SvrG): 

Together with e-folding time this leads to a 3-dimensional phase space. Insert x, y into 
the equations above to get: 

^ = ---+-il + 4^] (6.25a) 

'"^^ 'W-x^(l^4L\ (6.25b) 



dN 2 y dpm 2 V 9pm/ 
with constraint + = 1 for flat FRW space, the number of independent equations 
reduces to one: 

dx 3 3 f df \ , , 

Also, we know the universe if we know f{pm)- And for simple Cardassian terms, we 
can even get the exact and analyzable solution. 

6.2.5 Constraints for Cardassian &: Hyperbolic cosecant Cardassian 

Based on the general discussion above, we can make detailed analysis with certain 
form of Cardassian term. Yet the Cardassian term should not be arbitrary, as has 
been shown in the above algebra, and a potential candidate must obey at least three 
constraints IGD: 



(1) As for the baby universe with pcHtfl P, total energy density g{p) must reduces 
to p; or Cardassian term f{p) vanishes. 
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(2)Cardassian term is dominant recently with redshift z ~ 0(1), when gi^p) differs 
from p significantly, or /(p) differs from null significantly. 



(3)The speed of sound > must be positive. 

The constraints is in fact powerful. Only those that satisfies the first condition 
is worth further algebra during which the second and third condition proceeds to be 
tested. Three set of reasonable Cardassian terms have been introduced before: 

(l)Power-style Cardassian [53] 

/)' = /) + (6.27) 

where h and n is tuning parameters, and n < 2/3. 
(2)Polytropic and Modified Cardassian [55 j 



P — Pcard 



P — Pcard 



\ Pcard J 



(6.28a) 



1 + 



(Pcard \ 
Pm J 



9^1 ^ 



(6.28b) 



where q and v is tuning parameters; Pcard a constant density value for reference. 
(3)Exponential and Modified Cardassian [6 Ij 



p' = pexp(^^y (6.29a) 
p' = {p + Pcard) exp (^£^\ " (6.29b) 

\P + PcardJ 



where q and n is tuning parameters. 



Next we will discuss a new Cardassian term, namely hyperbolic cosecant Cardassian 
term: ^ 

p' = p(l + 6-csc/i(— ^ ) ) (6.30) 

V \ Pcard J J 

where b, n are tuning parameters and 

f[p)=ph.csch(-^\ (6.31) 

\ Pcard J 

Hyperbolic cosecant function differs very much from double exponential potential 
function. Before further discussion, it's a convenience to review the properties of hy- 
perbolic cosecant function. 
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6.3 Fundamentals of Hyperbolic Cosecant Function 

Hyperbolic sine sh{x) and hyperbolic cosine ch{x) are the two basic functions which 
are orthogonal to each other in the hyperbolic function space, and are defined as: 

p-^ p *^ p^ _|_ p ^ 

sKx) = ch{x) = (6.32) 

And definition of hyperbolic cosecant function is: 

1 2 

csch(x) = , , , = (6.33) 

^ ' sh{x) - e-^ ^ ^ 



Its one order derivative is: 

dcsch{x) 



= —csch{x)cth{x) (6.34) 



dx 

Where cth{x) is the hyperbolic cotangent function: 

eth(.) = = ^;±^ (6.35) 

sh{x) — e ^ 

The hyperbolic cotangent function can be expressed with hyperbolic cosecant function: 

cth{x) = Vl + csc/i2(x) (6.36) 

The inverse function of hyperbolic cosecant is denoted as Arcsch{x), and certainly 
csch{Arcsch(x)) = x. 

6.4 Csch Cardassian for unitary matter and radiation 

We treat matter and radiation together for the matter term, denoted 
the total energy density of hyperbolic cosecant (hereafter donated as Csch) Cardassian 
universe is, as above: 

p' = Pm(l + b- csch {y^^ (6-37) 

with the Cardassian term: 

f{p) = pmb ■ csch (-^) (6.38) 

\ Pcard / 

n and h are tuning parameters, which need observational constraints. Pcard is a constant 
density value for reference, as above in exponential Cardassian. when pm oo, p' — >■ 
Pm, thus the modified Friedmann equation returns to the normal form, and 

^ = h- csch [-^^ " -nh- csch f-^) " • dh (IH^ " . " (6.39) 

OPm \ Pcard/ \ Pcard J \ Pcard J \ Pcard J 
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Hence from (I6.18P we get the pressure that: 



df „ 

P = PmT^ / 



Pmb-csch 



Pcard 



-npmb-csch 



-hnpm ■ csch 



Pn 



Pcrit 
Pm 



Opm 

■cth 



Pn 



.(in 



cth 



Pcard J \ Pcard 
Pm \ I Pm 



-pmb-csch 



Pn 



^ Pcrit J \ Pcrit J \ Pcrit , 

the Friedmann quation and Raychaudhuri equation in Csch cosmos is: 



SvrG 



Pm + Pmh ■ csch 



( In 



\ Pcard 



Pcard 

(6.40) 
(6.41) 



and 



H = -47rGp^ 1 + 



dl_ 

dpm 



-AnGpr, 



1 + b • csch 



( 111 



\ Pcard 



nh ■ csch 



(^Ycthf 

V Pcard J V 



Pn 



Pcard 



f Pm V 
\ Pcard J 



(6.42) 

And the energy conservation equation p' = —3H{p' + p') is equal to the conservation of 
particles' number pm = —3Hpm , After the algebra above, we have known the kinetic and 
dynamical equations of Csch Cardassian cosmos when matter and radiation are treated 
together. Then we will research into its phase space and autonomous equations. We 
define: 



1 = ^^[Pm + fipm)] 
= ~l~ ^card 



(6.43) 



and introduce new variables = SvrG: 



K-VPn 



(6.44) 



thus ()6.39p can be rewritten as: 
dp 



1 + 



y 

bx"^ 



2 \ 2 



• Arcsch 



bx'^ 



(6.45) 



According to (j6.25p . we get the autonomous equations: 



dx 
dN 



3 3 / df 
2 2 V dpm 



3 3 
-x + - 
2 2 



n- 



1 + 



x^ 



.2 \ 2 



6x2 



• Arcsch 



6x2 



(6.46a) 
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dN 2 y dpm 2 V ^Pm 



- • 1/ 1 + ] ■ Arcsch- „ 



y 



, , y^ 

x^ 



1 + ( ^ ) • Arcsch^ 
ox^ / ox^ 



(6.46b) 



For flat FRW universe, rJm+f^card = 1) and there's the constraining equation x^+y^ = 1. 
So ()6.46p can be reduced into: 



dx 
dN 



3 3 
--X + - 
2 2 



1 1 , / /l-x2\2 1-x2 

- 1 - - 1) • Wl + I , . I -vlrcsc/i- 



3 3 / /I — x2\^ 1 — x^ 

--x^ — 2^(-'- ~ ' \ I ~ — ~ ' Arcsch- 



(6.47) 



bx^ J bx'^ 

This is the very difl'erential autonomous equation for Csch Cardassian cosmos with 
unitary matter and radiation. 



6.5 Csch Cardassian for separate matter and radiation 

6.5.1 Dynamical Equations 

In this section the energy density of matter and radiation will be treat separately, 
so the total energy density is: 



p' = {Pm + P-y) 
thus the Csch Cardassian term is: 



l + b.csch(P^^^X' 

\ Pcard J 



f{Pm, P-f) = f{Pm + /O7) = b{pm + P-y) ' CSch 



Pm + P-r 
Pcard 



(6.48) 



(6.49) 



As for Csch function, the symmetric property csch{pm, p^y) = csch{pm + Py) holds,so 
the matter term and radiation term are symmetric mathematically to each other. 

- - b-csch (p- + pA ^-nb-csch (p-^ + pA \cth f Pm + PyY f Pm + P.' 



dpm dp., 

\ Pcard 

According to 16.81 the total pressure is: 



Pcard 



Pcard 



Pcard 

(6.50) 



> Py ^ df df 
3 dpm 3 ' dpj 
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= y + (Pm + 3P7)^-/ 



-(Pm + tpj)'n'b ■ csch ( J . ci/i 

" \ Pcard J 



m 

Pm ~\~ P"i 
Pcard 



+ P7 \ Pm + Pj 



Pcard J \ Pcard 

Pm + P7 ' 
Pcard 



-(Pm + P^)b{pm + Pt) • CSC/I 

3 3 \ Pcard J 

- {pm + ■ nb ■ csch ( . cth ( ^i^y . ( ^^^y (6.51) 

3 \ Pcard J \ Pcard J \ Pcard J 

With Priedmann equation and Raychaudhuri equation being: 

^f' = ^[Pm + P7 + /(Pm,PT)] (6-52) 

and 

4 df 4 df 

H = -AiTG{pm + 0P7 + Pm-^ ^ qP7;h— ) (6-53) 

3 dpm 3 ' dp^ 

Here the energy conservation means the conservation of particles' number and the 

radiation conservation: 

Pm = -3-ffpm (6.54) 

p^ = -AHp^ (6.55) 

Now we have got the dynamical equations for Csch Cardassian cosmos with the radia- 
tion term and matter term treated separately. Then we'll investigate the autonomous 
equations in its phase space. 

6.5.2 Phase Space and Autonomous Equations 

Prom Priedmann equations, we define: 

1 = 3^[Pm+P7 + /(Pm,P7)] (6.56) 

= ~l~ ^7 ~l~ ^card 

where ^jn^^^,^card is the relative energy density for matter, radiation and Csch Car- 
dassian term respectively. Further more we introduce new variables with = SttG: 

x= — ;= — , y = — ;= — , z = —p^ — , iV=lna o.57 
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Thus ()6.50p reads: 



^'^ ^-Z^^-^- flTZ^^-^) -^^csc/i— 5^-^ (6.58) 



For flat FRW cosmos with 0^ + + ^^card, or + + = 1, combined with 
the symmetry, fipmiP-y) = f{Pm + P7), according to (|6.17p . we have the autonomous 
equations: 



+ ( + 2xy^ 



z2 / /l_^2_y2\2 l_^2_y2 

n—:: o\l+{-r7-^ ^ ■ Arcsch 



2;2 _j_ y2 W \^ ^(•3,2 _|_ y2^ J Ij^j,2 _|_ y2 



l-x2-y2 / /l_a;2_y2\2 l_3.2_y2 

n • ^ ^ — • 1 / 1 + ^ • Arcsch 



2-2 _|_ y2 W \^ ^^^2 _j_ y2^ J ^^2;2 + y2 



6.5.3 Critical Points 



From the above analysis, we get the autonomous equations in Csch Cardassian 
for separate treatment of matter and energy density. The following work is to get 
the critical points, and we will show there are three and only three critical points, 
namely:(0,0)(l,0)(0, 1). 

Obviously, 1 > x > 0,1 > y > and 1 > x^ + y^ hold, firstly let's introduce a new 
variable: 

^ ^ 1 - (x^ + y^) ^ Q 

X^ _j_ 

and formulate a new function of: 



= (j)- n</>v^l + 02/62 . Arcsch-j- (6.61) 



For (I6.6ip . when (x,y) — > (0,0), (j) — > +00, rupy^l + ^2/52 _^ Arcsch^ — > 0; and 
when (x,y) — > (1,0) or (0,1), i;^ ^ 0, n^y^l + ^2/52 g, ^rcsc/i| 00, The limit 
of ^ is impossible to be read directly in both approximation cases. So we transform ^ 
and employ the familiar Los'pital theorem (prime operator (')refers to derivative of the 
dependent variable to the independent variable): 

Arcscht 
^ = 6 , ^— 
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Arcsch 



= ct>- 



2J.2 



^ b 



(6.62) 



Based on this equation, let's consider two approximation cases: 
(1) When (x,y) ^ (0,0), 4> ^ oo: 



+ 



+ 



(6.63) 



nb ■ (j) 



So the infinitesimal Arcsch^ is magnified by the infinity n(f)^l + cjp' /b"^ to get an unitary 
infinity n(f)^l + (f>'^ /b'^- Arcsch^, which is coincidently of the same order with the infinity 
(f), who will conduct a minus operation with it. The minus of two infinity of the same 
order can leads to conservance ! And we put (f) back to the original equation, and when 
{x,y) — (0,0), we have: 



dx 
dN 

> 



-x^ + 2xy^ 



-x-^ + -xy^ 



• (1 - nb) 



1 - (x^ + y"^ 
x^ + y^ 



(6.64) 



(1 - nb) 



1 - (x2 + y2 



_|_ y2 



^x-(l-|n)[l-(x2 + j/2)] 







It's notable that ^ is conserving itself if we put 1 — n6 = 0. However, this constraint is 
too strong and we will have much less room left for parameter tuning. Considering there 
is an infinitesimal factor waiting to multiply we'd better realize the conservation in 
the next step, rather than set 1 — n6 = in a hurry in advance. Hence (0, 0) is a critical 
point of this autonomous system. 



(2)When {x,y) (1,0) or (0, 1), ^ 0, in 



+ n 



2j>2 



the numerator is higher order infinitesimal compared to the denominator, so the limit 
is zero, and: 



+ n 



2j,2 



~S3 I f 



+ 



ncp 
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whose limit is also zero. Now let's go back to analyze the original equation: 



dy 
dN 



2y + 



-x'^y + 2y^ 



3 

2y + -x'^y + 2y^ ~ 

Hence (1,0),(0, 1) are both critical points of the autonomous system. 



(6.65) 



All in all, (0,0), (1,0), (0,1) are three critical points of this autonomous system. 
And because the critical points are numerically achieved via solving the group of two 
three-order algebraic equations, according to Gauss's basic theorem of arithmetics, we 
could have only three roots in the complex domain. So the three critical points achieved 
via approximation analyzing are the only ones. 

6.5.4 Stability Matrix for Autonomous System 

Now we'll make phenomenological analysis of the stability of the critical points. 
Set X = ■^,Y = and the perturbation matrix will be: 



M 



dX dX 

dx dy 

dY dY 

dx dy 



(6.66) 



The four elements of this Jacobian matrix is: 
dX 3/9 2 



1 + 



1 - (x^ + y 

^.2 _|_ y2 



2\1 



-x^ + 2y^ 



n- 



+ ( +2xy Jr^-^^2 + y2)2 



1 - (x^ + 
x^ + y2 



2x 



'1 + 



1 — (x^ + y2) 
6(x2 + y2 



2^\2 



Arcsch 



Arcsch 



,'3 3 ^ o\ (l-(x2+y2))2 

+n [ -x^ + 2xy^ ' ^ ^ ^ " 



b{x^ + y^) 

2:r 



Arcsch 



1 - (x^ + y^ 
6(x2 + y2) 

l-(x^ + j/^) 
6(x2 + 2/2) 

1 - (x^ + y2) 



i^'+y')' &Vl + ^(i4p-Ij^ 6(x2 + y2) 



2x 



(,2 + ,2)2 + ^-y') + + 2xy^) "(^^1^ ^'-''^ 



and 



dX 
dy 



4xy 



1 + 



1 - (x^ + y 
x2 + y2 



2\-\ 



x2 + y2 



6(x2 + y2) 



6(x2 + y2) 
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1+ 



6(x^ + y^) 



Arcsch- 



+ -x"* + 2xy^ n 



(l-(x2+y2))2 



2y 



I - {x^ + y^^ 

2 I ,,2^ 



Arcsch 



2y 



x"" + 2a;y^ + -a;'' + 2xy^ n 



2y 



and 



9y 



1 + 



2~ ■ — - y --(^2 + ^2)2 



(6.68) 



^2 + 2/2 



1 - (x2 + y2^ 
a;2 + y2 



-3xy ■ n 



1 + 



2 _L o,2\ \ 2 



1 - (x2 + y' 
6(a;2 + y2) 



1 - (x2 + y2^ 
6(^2 + y2^ 



Arcsch 



+ l2^2/ + 2y j" (^2+^2)2 



9 -; • 



1 + 



6(^2 + y2) 



Arcsch- 



6(x2 + y2) 



,'3 2 3\ (l-(x2+y2))2 



2a; 



Arcsch 



l-(x2+y2) 



(^' + ^')' 6(x2 + y2) 



6^ ^ x'^-\-y'^ 

/So o s\ 2 ..A 2,T 



(6.69) 



and 



ay / 3 2 2 



1 + 



i-(.r^ + r) 



x2 + y2 



l-(x2+y2) /i_(^2 + y: 



1 + 



r.2 _L „2\\ 2 



6(^2 + y2) 



Arcsch 



6(^2 + y2) 



+ (2xy + 2y Jn^^2^^2)2Yl+^ 6(x2 + 2/2) ) Arcsch ^^^2 + ^2) 



+ ( ^x^y + 2y^]n 



(l-(x2 + y2))2 



2x 



1 - (x^ + 
2 I „2^ 



Arcsch 



2y /3 



(x2 + y2)2 



x^y + 2y3 ) + I -x^y + 2^^)71 



2„, , r,„.3 



2y 



(x2 + y2p 



(6.70) 
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Chapter 7 

Virial Collapse of Power-style 
Cardassian 



7.1 Applications of Virialization in Cosmology 

The spherical collapse mechanism, which was put forward by Gunn and Gott [82] 
|83j . is an effective tool in understanding the consequence of the density inhomogene- 
ity. An overdensed volume will slow down from the co-expansion, until the expanding 
velocity vanishes. Then it starts to collapse over rotating, and potential energy trans- 
lates into rotating kinetic energy. The special state when the expanding just stops and 
collapse is to start is dubbed turn around state. Ultimately the collapse ends up with 
stable virial structure. 



We can apply the virial theorem to the above process in the light of energy conser- 
vation. A general statement of virial theorem says that, during a long period of time, 
the opposite value of the average kinetic energy of a particle group equals to the virial 
force acting on it: 



T -- 

2 



where T is the total rotating kinetic energy, Fi and fi refer to the external forces and 
location vector of the ith. particle. For conserving system, 

f = \Y,(y~^)-n (7.2) 

i 

where V is the potential of the particle group, V is the vectoring gradient operator: 

Vi = ^i + ^i + ^A; (7.3) 

dxi dui dzi 
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The total potential energy, the total virial potential energy and the static radius 
of the virialized state at the turn around (when the expansion stops, and collapse and 
rotating takes place) are denoted as Uta, Umr and Rmr respectively. Supposing the 
system's potential energy takes the analytic form of U, then the kinetic energy after 
virialization, which depends on U, is: 

The energy conserves, so after the virial collapse(T=0 when collapse starts): 

RdU 



^+2 9i? 



Uta (7.5) 



The analytic expression of the potential energy term U can be deduced as follows: 
Supposing pq is the energy density for the virial component, is the total potential 
acts on a certain point in the sphere, and the total potential energy will be 

^ = 11 PQ'^QdV (7.6) 



If 4 o 



In a homogeneous sphere, pq satisfies: 



r 



Pq + 3{1 + wq)-pq = (7.7) 



And the total potential in the sphere for a certain point is: 

^2- 



<l>Qir) = -2^G(1 + 3wq)pq (^R^ - (7.^ 



Where R is the radius of the sphere, r is the distance between a point and the spherical 
center. 

7.2 Virialization Power-style Cardassian 

Now let's consider the virialization of power-style Cardassian [53]. Supposing that 
Cardassian term and matter term both participate the virialization process, i.e. 

p = p^ + bpl„ n<2/3 (7.9) 
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7.2.1 Virialization of Matter Term 

EOS parameter for matter term is if = 0, and the distribution of potential in the 
sphere is: 

$(r) = -27rG(l + 3 x 0)/9„ ( - ^) (7.10) 



The total potential energy in the sphere of radius R is: 



„2 



U^ = - I ■ 2TrG ■ 2TTr^ ■ [ R^ - — ] dr (7.11) 







16 2^ 2 7-,5 

The rotating kinetic energy associating Um is: 

= f ^ = -l^'GplR' (7.12) 

Obviously Tm < in (I7.12p . Kinetic energy cannot be negative, so the virial collapse 
of matter term alone is forbidden. 

7.2.2 Virialization GF Fluid/Cardassian Term 

Now let's have a look at the behaviors of Cardassian energy. Power-style Cardassian 
is perfect fluid with EOS parameter w = n — 1 and the distribution of potential in a 
virializing sphere is: 

$(r) = -27rG(3n - 2)(6p;^) (i?^ _ (7,13) 
So the total potential energy in the sphere is: 

U =^ j p$(r) • Anr^dr (7.14) 

= -l£pm- 2vrG(3n - 2)ibp^) (^R^ - • ATrr^dr 

= ^G7T'R'{2-3n)bp'^ 
and the associated virial rotating energy is: 

n^r = ^^ = ^Gn'RU^ - in){p^,,.r + hpl,^,„f = hj^ir (7.15) 
Energy is conserving during the virialization process. 

! X ''±Gi?Rl„{2 - 3n)(6p^V) = ^G7r2i?5j2 - 3n)(6p^",,J (7.16) 
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reduced to be: 



The mass does't lose, so 



Put mHl int dTlTl to get: 



7 



Pm,vir / ^ta 



Pm,ta \R 



(7.18) 



Rta J V Pm,ta J \ Rta J \ Rvir J \ Rta J 7 

Hence the ratio of the stable virial radius with turn around radius is: 

= ^-^"/? (7.20) 

Rta V 7 

The upper calculations constrains parameter n for twice, namely kinetic energy being 
positive and collapse radius becoming smaller than initial: 

T>0 ~2-3n>0 (7.21a) 



Rta V 7 

So, after virialization GF fluids collapse to a steady sphere with certain radius which 
depends on the initial radius and the Cardassian parameter n. And no singularity is 
generated. 

7.2.3 Unitary Virialization of Matter Term and GF Fluid Term 

Although matter term pm cannot virialize on itself, howerever, as it is shown below, 
due to the excellent virialization property of GF fluid, the matter term pm can virialize 
together with GF fluid Pcard- This is pretty similar to the fact that high energy 7 pho- 
ton cannot decay into electron-positron pairs when isolated, but will decay if located in 
Column field nearby the nucleus. 

As analyzed in chapter 6, that the total pressure of pm and pcard is p = (n — l)6p^, 
so the EOS parameter is: 

in — 1)1)0^ 

w = 7-^ 7.22 

Pm + bp"^ 

The distribution of potential generated by pm and Pcard in the virial sphere is: 

Hr) = -2nG (l + ^^^^^|^) (P^ + hp^) (r' - (7.23) 
-2TiG{pm + Kn-2)bp%) 
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Pm and pcard leads to the potential: 

U=l ri-27rG){pm + 3in-2)bp^)dr 

^ Jm 



= ^'K-'GR\p^ + bp^J{{2 - 3n)bp^ - p„) 
the associating rotating energy: 

T = ^TT^GRHpm + bp^,){{2 - 3n)bp^ - p^) = 

Thus, T > is possible if b, n are properly tuned. Due to energy conserving 
T 16 

2 ^ -^^"^^^liriPmMr + Vm,mr)((2 " 'in)bp'^ - Pm) 

= ^^-n^GRUPmM + bp^^ J{i2 - 3n)bp^^ - p^) 
The mass is conserving 

Pm,vir / ^ta 



Pm,ta \ Rv 

Combine it with energy conserving equation: 

7 / Rvir \ ^ {Pm,ta + bpl,^ta) ((2 " 3n)6/9;5, - Pr, 



2\Rta J {Pm,vir + ^p^^^i^) ((2 - 3n)5/9^ - Pr, 

_ (2 - 3n)bV^^,, + (1 - ?,n)hp^^^l - p^,,, 
(2 - 3n)6V^"„. + (1 - 3n)6pr-l - pI 



m,vir rm,vir 

Set the ratio = x, and via the mass conservation 



7 5 (2 - 3n)62p^«^^^x6« + (1 - 3n)bpl%x^-+^ 
-X — 



,2 . 
rriMr 



2 (2 - 3n)6V^%,, + (1 - Mbp':;^X " 

Make deduction of x^ to (fr29]) to get 

^ {(2 - 3n)6Vm,«r + (1 - 3n)6p:^+V " P^vir] 
= (2 - 3n)6V^^..^'"-' + (1 - 3n)6/,;^i,x3"-2 - p; 
Introduce a new variable g = Es^s^^^ = ^ and 

Pm.vir Pm,vir 



>,2p2n ^ 2 2 



rm,vir ^rm,vir 

Insert (17:31]) into (f7:30ll to get 

7 r '- - - _2 2 



- {(2 - 3n)q p^ ,,^,. + (1 - 3n)gp„,„i,. - p 



m,vtr 
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= (2 - 3n)q'pi,^,i,x^'^-^ + (1 - 3n)gp^,,,,x^"-^ - pt^„,x (7.32) 

Obviously Pmvir > holds, and deduct it from both sides of the equation (|7.32p . and 
we get the relation of x and q: 



-{(2-3n)g2 + (l_3n)g-l} 

= (2 - 3n)g2^6"-5 + (1 - 3n)gx3"-2 - x (7.33) 
We can also get the relation z^r and x via ()7.33p . Currently, Friedmann equation reads: 

SttG 

Hi = ^{Pm,0 + bpln) (7.34) 



I.e 2 

Pcrit,o = = Pm,o + bpl,fl (7.35) 

i.e 

Pm,0 Pm,0 

i.e 

-^-l = q (7.37) 
insert it into x — q relation (17.330 and one will get the q — Zyir relation. 



7.3 Conclusions and Evaluation 

New work gathers in Chapter 6 and Chapter 7. Firstly we make a comprehensive 
analysis of the existing work, pay particular attention to Gondolo and Freese's idea of 
treating Cardassian energy term as relativistic perfect fluid(GF fluid), point out that 
a potential Cardassian term should meet three conditions, and review three existing 
Cardassian terms (power style, polytropic style and its modification, exponential style 
and its modification) , and eventually put forward the newly found hyperbolic cosecant 
Cardassian. 



Then the Cardassian dynamical equations are introduced generally and logically 
under GF fluid scenario, together with the flowing process of constructing phase space 
and differential dynamical systems from Friedmann equation. Hyperbolic cosecant Car- 
dassian term is employed for concrete computing. The analysis proceeds in two cases, 
namely a unified description of matter and radiation energy density (case 1) and a sep- 
arate description of matter and radiation terms (case 2). Formalism of case 2 is more 
exact at the expense of more complicatedness, and due to the mathematical symmetry 
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of matter term and radiation term in hyperbolic cosecant function, the differential dy- 
namical equations are considerably simplified. Phase space and dynamical systems for 
both cases are achieved. When we calculate the critical points for case 2, amazingly 
interesting behaviors of self-consistency and auto-normalization are exhibited, which is 
a strong support for the new model, along with a forever positive sound speed. 

The process of virial collapse in Cardassian cosmos is analyzed. Power-style Car- 
dassian term is employed for its simplicity. Calculation declares that virial collapse of 
matter alone is forbidden. Yet Cardassian has excellent ability for virial collapse, after 
the virial collapse ending up with a stable sphere, the ratio of the ultimate radius to 
the original radius depends on the adjustable parameters in Cardassian term. And, the 
mixture of GF fluid and matter could conduct virial collapse, the ratio of the ultimate 
radius to the original radius depends on the adjustable parameters in Cardassian term, 
too. 

The creative work in this thesis incudes the introduction of hyperbolic cosecant 
Cardassian, which is the fourth ever available Cardassian term ever found, after the in- 
troduction of exponential style in 2005. This work, along with the calculating of power 
Cardassian virialization, enriches the research of Cardassian cosmology. 
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